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In this paper, we discuss the effect of the arbitrary spin impurities to the 

spin-1/2 and spin-1 XXZ model. The effect of ground state, the free energy, 

the magnetic susceptibility, the specific heat and the Kondo temperature are 

given by using the thermodynamic Bethe ansatz equation. 
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I. INTRODUCTION 

Behavior of a quantum chain will be quite different when impurities embedding in the 
system. With the method of bosonization and renormalization techniques, Kane and Fisher 
proved this result when they studied the properties of a potential scattering center in a 
Luttinger liquid [0]. After their work, many papers have paid their attention to the problems 
of local perturbation to a Luttinger liquid, especially the Kondo problem in the system [2-5]. 
It is well known that the spin dynamics of the Kondo problem is equivalent to the dynamics 
of the spin chain with magnetic impurities [^]. 

The Heisenberg chain is an important model in the study of integrable system. Firstly, 
Andrei and Johannesson considered a spin-s (s > |) impurity embedding in the spin-| 
Heisenberg chain with periodic boundary condition 0. Then, Lee and Schlottmann gener- 
alized the problem to arbitrary spin ||^ . However, it is found that the open boundary theory 
is better to deal with quantum chains with impurities, since the periodic boundary condi- 
tion would lead to the presence of some unphysical terms in the Hamiltonian to maintain 
the integrability 0,|l3l- In resent research, in the frame of the quantum inverse scattering 



method , Wang studied the exact solution of the spin-| XXX open Heisenberg chain with 
two arbitrary spin impurities 0. Subsequently, Hu and Pu considered the spin-i XXZ chain 
with two spin-| impurities |]TU|. In the previous paper [|TT|, we obtained the eigenvalues of 



the Hamiltonian and the Bethe ansatz equations for both the spin-| and spin-1 XXZ chains 
with two arbitrary spin impurities. 

The low-temperature thermodynamics was first proposed by Yang and Yang. In their 
series papers, they studied systematically the ground state, the magnetic susceptibility and 
the specific heat of the spin-| Heisenberg chain [p!^-[l^. Directly, this method was general- 
ized to other models, such as the Hubbard model etc [|15|-|1^. In the Ref. |T8|-p2|, Babujian, 
Kirillov and Reshetikhinthe and others studied thermodynamics of the arbitrary spin XXX 
and XXZ Heisenberg chain respectively. In the present paper, we study thermal properties 
of the spin-i and spin-1 XXZ chains with two arbitrary spin impurities by using the standard 



method at low temperature. 



II. THE MODEL 



The integrabihty of the spin-^ and spin-1 XXZ Heisenberg chains with two arbitrary spin 



impurities have been discussed in the Ref . ||1T| . Hamiltonians of them can be represented by 

(1) 



—i d t(u) 

H = r--- ; ^ : ^ -— — X 



M=0 



2p^(0)pd(c„)prf(cb)tr,i^+(0) du 

where t{u) is the transfer matrix, p(0) and pd{ci) {i = a, b) are constants appearing in 
unitarity relations of the i?-matrix in the bulk and the boundary respectively, K^ is the 
dual reflection matrix, and Cj are coupling constants. 

Applying the Hamiltonians H on the multi-particle-state Bethe vector \Q), we obtain 
the eigenvalues as follows 

E = Eb + E, + E^, (2) 

where Ef, , Ei and Eoo stand for the energy of boundary , impurities and the bulk respectively. 

To the spin-i chain , we have 

^ ^ i sinh(^+ - ^ i 

2sinh^+sinh^ 2sinh(2?7) 

^» = ~9 H H coth(rCfc + SiT] + -) , 

■^ r=±l k=a,b ^ 

^ ^ _^i_coshr/ _ ^ 2isinhr] 

°° sinhr] ~[ cosh(2fj) — cosh?] 

where Sj is the arbitrary spin of impurities; to spin-1 chain case, we have 

zsinh(2?7) icoshrj 

^ ^ ~2sinh(e+ + f) sinh(e+ - f ) ^ sinh(3r/) ' 



„ v^ v^ ^ sinh(2rcA; + 2sir] + t]) 



T^i fct^b 2 sinh(rcfc + s^r/) sinh(rcfc + SiT] + r]) 



^ ^ ^^, i sinh(3r/) _ ^^icoshj^ _ ^ 2i sinh(2r7) 

°° sinh?7sinh(2?7) sinh(2?7) ~^ cosh(2fj) — cosh(2?7)' 



where ^^, ^ and t] are some free parameters and Vj satisfies the Bethe ansatz equation 
sinh(t>j — ^+ — |) sinh(t>j + ^ — |) sinh(t>j + |) cosh(t>j + |) sinh^^(t;j + srj) 



sinh(t>j + ^+ + |) sinh(t>j — ^ + |) sinh(t>j — |) cosh(i;j — |) sinh {vj — srj) 



sinh(f j + rck + Sii]) 



^ sinh(t>j — Vi + 7]) sinh(t>j + Vi + r]) 



riii kiab sinh(t;j + rCk - Sif]) fj^ smh{vj - Vi - rj) smh.{vj + Vi - rj) 



(5) 



where j = 1, 2, 



M with M being the number of down spins. 



With the Bethe vector \il), according the Ref. [T^JT^, we can also obtain the eigenvalue 
of the total spin z-projection 



Ns + 2si- M 



(6) 



III. THE GROUND STATE 



To obtain precise forms of the eigenvalues of Hamiltonians, denoting 



^^{vj) = 2 tan ^ 



7 
cot — tanhf,- 

2 \ 



and taking logarithm of the Bethe ansatz equation (5), we have 



(7) 



M 

+ E E "^^sA^J + rCk) = 27rlj + J2 {^jivj - v^) + ^^{vj + Vi)} , (8) 

k=a,br=±l i=l 

where the transformations t] ^ i'y, C, ^ iC, and C,~^ -^ iC,~^ have been used, and Ij is an 
integer. Taking the thermodynamic limit ( A^ ^ oo ), the distribution function of particles 
can be considered as continuous one. So the Bethe ansatz equation can be written as 



ZTT 



*L7(^) + ^d'i^) 



1 f^ 

^y_^rf'^'^27(^-^)p(^) 



(9) 



where the distribution function is defined by ||T2 



p{v) 



dij (v) 
dv 



(10) 



and 

k=a,br=±l 

Here the 6{v) stands for the effects of open boundary condition and magnetic impurities. 
Applying the Fourier transformation on eq.(9), we have 

sinh(7ru; — 2s'yuj) 1 j cosh(|u; — 'ju) sinh(|c<j~) 



+ 



2 sinh(7rco' — 70;) cosh(7c<j) 2N \ sinh(7rto' — 70;) cosh(7c<j) 
cosh[(7r - ^ - ^+)uj] sinh[(^+ - ^ + -f)u] 
sinh(7rci; — 70;) cosh(7c<j) 

2 sinh(7ru; — 'yuj) cosh(7u;) J ' 
where the B{uj) = 4cos[ca + Cb)uj] cos[(ca — Cb)uj]. 
The distribution function p{v) can be obtained by 

/oo 
p{uj)e-^'''^duj. (12) 

-00 

The ground state energy per site for our systems can be written as 



E^E^E^E^^ (13) 

N N N N ^ ^ 

(i). when s = |, after some calculation, we have the ground state energy of the bulk term 

^ , cosh 7 ^^ f°° , sin(7) . . 

i^oo = -N—r-^ + N / dv ^^^ 14 

smh7 J -00 cosh(2s7j — cos(2f J 

It is exactly same as one in the periodic case |T^ . The contributions of the boundary ( Eb 

) and the impurities ( Ei ) are 

^^ ^ sinh(e+-0 ^ 1 

2sinh^+sinh^ 2sinh(2r7) 

^^ sinh[(e + - e + 7)cu] cosh[(7r - ^+ - 0^^] 



sinh(7rti;) cosh(7Co') 
cosh(7rcij/2 — juj) sinh(7r(x;/2) , , 

sinh(7rco') cosh (70;) 



-00 
^i= n Yl Yl coth(rCfc + isa + — ) 

^r=±lfc=a,fe ^ 



°° 4 sinh[(7r - 2si'y)uj] cos[{ca + Cb)uj] cos[(ca - Cb)uj] ^ , . 

-00 sinh(7rco') cosh(7Ci;) 
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(ii). when s = 1, similarly, we have 

E = N ^^^^^^'^^ _ 2ivS^^!i^^M 
"" sinh7sinh(27) sinh(27) 

,^ r°° , sinhfTTcu — 27Ct;) 



■'oo 



2 sinh(7rci;) cosh(7Co') ' 
sinh(27) 2 cosh 7 

'' " ~2sinh(e+ + f)sinh(e+-f) ^ sinh(37) 

^^sinh[(e+ -i + 7)cu] cosh[(7r - e+ - 0^] 



sinh(7ra;) cosh(7Co') 
cosh(7rci;/2 — 70;) sinh(7ru;/2) , , 

^^ ^T7 ^ TT — ^ (^^'^ 

smh(7rcc;j cosh(7ti;j 

_ v^ v^ isinh(2rcfc + 2isj7 + ^7) 

rt±i kt^,b 2 sinh(rcfc + isa) sinh(rcfc + isa + i-f) 

4sinh[(7r - 2sa)uj] cos[(ca + Cb)uj\ cos[(ca - Cb)uj] 

sinh(7rci;) cosh(7Ci;) 



IV. TBA AND PHYSICAL PROPERTIES 

It is well known that that the general solution of Bethe ansatz equation (5) in the 
thermodynamic limit lies in the complex plane and forms strings with the length n 0,^: 

^3 = K,j = K + ^\{n + 1 - 2j) ( J = 1, 2, • • • , n ), (19) 

where n = l, 2, ■■■, 00, and a = 1, 2, ■ ■ ■ , M^"'\ Here the number of n-strings M^") 
satisfies X^^i ^^'•"^ = M. Substituting the above string solution into eq.(5), we have 

2iV^:(A^) + ^„(A^) = 2vr/:+ ^ {y,_(A^ _ A™) + y,_(A^ + A^)} . (20) 

m,/3>0 

Here, 



c(a::) = T"'^"^^ ' ' 



+ H '^7("+2..+l-2i)(A^ + rck) 
r,k 

V5mn(A^ ± A™) = (1 — (5mn)'^7|m-n|(AQ ± AJ^) + 2$^(|m-n|+2) (AJ^ ± AJ^) 
H h 2$^(|„+„|_2)(A;J ± A)^) + $^(r„_„)(A;j ± X^) 
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Under the thermodynamic hmit, according to [|I4[, we obtain the following equation 



1 1 fOO 

Pn(A) = ^[«)'(A) + T^e'nW] - E / ^M..n^(A - /x)p„(/x), (21) 

27r 2N „>i^-oo 

where pn and p^ are the distribution functions of the particle and the hole respectively, and 

zvr 

The energy of per cite can be written as 

K 1 c°° 

Tf = TfE}-Y. rfA(V^:)(A)p„(A), (22) 

where 

{A^ cosh 1] "5 = 1 

sinh 
The entropy of per cite can be obtained after we analysis the distributions of particles and 

holes. It reads 
S _ 

n>l' 

(23) 



- = E / d\ {(p„(A) + p^::(A)) ln(p„(A) + p::(A)) - p„(A) lnp„(A) - p::(A) lnp::(A)} . 



The eigenvalue of the spin z-projection can be obtained with the so-called Bethe vector. 
With the distribution function, we have 



2 r°° 

^+M^i-T.^ Pn(A)rfA (24) 

iV ^^^ J-co 



s' 2 



In order to find the equilibrium distribution functions at the finite temperature T, we 
must minimize the free energy 

F = E-TS- HS^ (25) 

with 5F = subject to the constrain eq. (22). Then we obtain the thermodjTiamic Bethe 
ansatz (TBA) equation 

ln(l + ^) = ^{Hn - Im'i^)) + E ^nmHl + ^), (26) 



where symbols with a hat are results after Fourier transformation on them, and 

-j 2 coth{'yu) sinh[rnin{m,,n)'yuj] smh[nu — rnax{'m,n)'~-fuj] , , 

Anm[^) = r^-7 ^ , (27) 

^nm\^) — ^nm — TT T [Sn.m+l + Sn.m-l) ■ (28) 

2 cosh(7C(j) 
Define 

5„(cu)=Tln^. (29) 

Then, the above TBA equation can be rewritten as 

ln(l + e-(-)/^) = ^Hn - i(V^:)'(c.)) + Y. A.„.ln(l + e'-^^)/^) (30) 

Solving the in{^) from the above equation, we have 

ii{uj) = Tpln{l + e=""+i(^)/^) (31) 

in{uj) = Tj51n[(l + e^~"+^('^)/^)(l + e^'-^^")/^)] - npSn,2s (32) 



where p = ^ — h — r- From eq. (28), it is easy to find that 



2 cosh(7Cij) 



e. 



lim — = H. (33) 



n 



Taking into account (26) and (27), the free energy per site at equilibrium state is given by 

F 1 1 
_ = —Ef - 2tan(2n7) - H(s H Si) 



/oo 
rfA(V^:)(A)ln(l + e- 



27r ., 

n>l 

We can now obtain the specific heat in the condition T ^ 0, H = 0, and the magnetic 
susceptibility at the small magnetic field by using the following formulas 



dT \N 



-T 



d' (F 



T 
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(35) 

H 

(36) 



Substituting eq.(34) into eq.(35) and (36), we find that the former three terms are zero 
and the fourth term is exactly same as the periodic case which has been discussed in Ref. 



T5|,p!8|-p0| . Then, we shall discuss the last term of the eq. (34). As before, the 6'„(A) can be 



divided into the boundary and the impurity terms 






(37) 
(38) 



r,fc 



Then, the last term of eq. (34) can be rewritten as 
F 1 T 



N 



^,,^/y^&^mix)Hi+e-^^^y-) 



E 



'(V^^)'(A) 



1 T 

'N2^ 



E 

n>l 



- (^^)'(A) 



A:^«)(A)ln(l + e--W/^) 



1 T 

AT O^ ^-^ 



N2tt 



n>l 



1 T 

--— E 



°^ rf^l$lM(^^)'(^) in(i + e-"H/^) 
" t/u;l$lM(^^)'(^) in(l + e--~"H/^) 

'- {rn)'{^) 



(39) 



In this paper, we mainly discuss the effect of the impurities so here we only give the result 
about the last term of the eq.(39). For both spin-i and spin-1 cases. 



(^;)'(^) 



^ „, , sinh(2si7u;) sinh(7ru; - n7u;) 
2'nB{uj) r-TT r— tt ^ n > 2si 

, , sinh(n7Ci;) sinhfvru; — 2sj'yuj) 
27rB uj) ^^-4 — A^^ — r^-^ n < 2s,- 



(40) 



sinh('7rLi;) sinh(7Li;) 
where B^u) = 4cos[(ca + Cb)u!] cos[(ca — Cb)uj]. And 

sinh(7rc<j — w-fuj) 



(^:)'(^) = 2n 



sinh(7rtt;) 
2 cosh(7u;) sinh(7rc(j — w-fuj) 



2 ' 

s = l) 



(41) 
(42) 



sinh(7ru;) 

Substituting the results of (^^)'(co'), (ip^y^uj) into eq.(39), we can obtain the effect of the 
impurities to free energy and other physical quantities. However, the calculation is much 
more complicated than we expected. Fortunately, comparing with the periodic case, we 
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can obtain the low-temperature specific heat and the magnetic susceptibihty contributed by 
impurities as 



C 



2C? 



1/2 



1/2 



N 



X cosh 



2C° 
C{ = -ZT- X cosh 



Xl/2 



x\ 



N 

2x?/2 

TV 

2x1 



X cosh 



N 



X cosh 



{cg + Cb)7r 

27 

Cffl + Cb)7r 

27 

(Cg + Cb)7r 

27 

Cg + Cb)7r 

27 



cosh 



Cb)7r 



cosh 



cosh 



cosh 



27 

(cg - Cb)7r 

27 

(cg - Cb)n 

27 

Cg - Cb)7r 

27 



sm 



vrsi 



(2K+l)7r^ 
27 



cosi^j 



sm 



vrs,; - 



(2E"+l)7r^ 

27 



smlt^J 



sm 



TTS,; - 



27 



sm 



cos 

TTSi - 



27^ 

(2i^+l)7r^ 

27 



sinlf:^J 



(43) 
(44) 
(45) 
(46) 



therefor, the Kondo temperature can be given by 



rpk 

-^1/2 



cosh 



-1 



Tf = cosh"^ 




cos 



27' 



2 sin 



ITS 



{2K+l)Tr^ 
27 



2 sin 



TTSi 



{2K+l)n^ 
27 



(47) 
(48) 



where C° and Xs ^^^ the specific heat and magnetic susceptibihty of the infinity size system 

determined by the relation < 2sj7 — 2K7i < 2'k. One can 



0^,111, and ir=0, 1, 2, 



check the results eq.(43)-(46) to be sound based on the fact the Cj and Xi proportion to 
the shift of the state density at Fermi surface which is proportional to piiy) j pooiy) ( see 
eq.(ll)and (12) ), where Piiy) and Pooiy) are the density of the impurity and the bulk. If we 
choose Cg = Cb and Sj = 1/2, the results of eq.(43), (45) and (47) can recover those in Ref. 
T0| . In eq.(47) and (48), when c -^ ic, the Kondo temperature satisfy the relation T^ ~ 



cos ^[(cg + Cf,)7r/27] cos ^[{cg — 0^)71/27]. The result shows a crossover from an exponential 
law to a power law as the phenomenon in Ref. 0. 

V. DISCUSSION 

In this paper, we obtain the contribution of the boundary arbitrary spin impurities to the 
susceptibility, the specific heat and the Kondo temperature by using the standard method at 
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low temperature. The model that XXZ chain coupled with boundary impurities is relevant 
to the Kondo problem in a Luttinger liquid. And the effect of the impurities to the bulk 
can be regarded as a perturbation to the bulk. From eq.(43)-(46), we can see that the 
contribution of impurities strongly depend on the coupling constants Ci {i = a ,b) and the 
spin of the impurities Sj. Also from these equations, the Wilson ratio can be obtained as 
{Cl/Cs)/iCs/^s) = 1 ill this model. This fact can prove that the present model satisfies 
the theory of the Luttinger liquid. Besides the properties discussed in this paper, other 
properties of the present model are also worth investigating, such as the critical properties 
of the present model. And we shall study them in our further papers. 
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